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Abstract — Caire, Taricco and Biglieri presented 
a detailed analysis of bit interleaved coded modulation, 
a simple and popular technique used to improve system 
performance, especially in the context of fading channels. 
They derived an upper bound to the probability of error, 
called the expurgated bound. In this correspondence, the 
proof of the expurgated bound is shown to be flawed. A 
new upper bound is also derived. It is not known whether 
the original expurgated bound is valid for the important 
special case of square QAM with Gray labeling, but the 
new bound is very close to, and slightly tighter than, the 
original bound for a numerical example. 

Index Terms — Bit interleaved coded modulation, 
BICM, expurgated bound, probability of error 



I. Introduction 

A comprehensive study of BICM is presented 
in [1]. There, in addition to an information theoretic 
analysis of BICM, a detailed analysis of the probability 
of error is presented. In the error analysis of BICM in 
[1], various upper bounds and approximations to the 
probability of error are derived, notable among which 
is the expurgated bound. In the first half of this paper, 
counter examples are given for the two theorems in 
[1] leading to the expurgated bound. Consequently, the 
validity of the expurgated bound in [1] is questionable. 
The second half of this paper focuses on the important 
and practical case of square QAM constellations with 
Gray labeling. For such cases, an alternate upper bound 
is presented. Numerical results are given for 16-QAM 
and 64-QAM and a rate-i convolutional code. For these 
examples, the new bound is nearly equal to, and slightly 
tighter than, the expurgated bound of [1]. The reader is 
referred to [1, Sections 2 and 4] for notation. 

II. TWO COUNTER EXAMPLES 

Counter example to [1, Theorem 1]: Consider 
the constellation in Figure \l\ It is similar to QPSK, 
except the point which is on the unit circle but 
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Fig. 1. 4PSK - The shaded regions correspond to the region not 
covered by the expurgated union bound. 



closer to one of its neighbors. The labeling // is chosen 
to be Gray labeling. For simplicity, assume d = 1. 
For concreteness, assume that the message to be sent is 
6 = 0, and that it is to be sent over the first label position, 
i.e. S = 1. Further, let U be so that the signal labels are 
not complemented. Then, an element of Xq = {00,01} 
is transmitted. Suppose x = 00 is transmitted. Let z 



> 



11. As in [1], let r M = {y : ||y-X||5 
Clearly r x z C T x z( i) U T x , (2 ) where z 
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10. So, by [1, Theorem 1], T x z can be neglected 
in the union bound. This leaves only T x B (aj in the union 
bound, but that fails to cover some part of the pairwise 
error region. In particular, referring to Figure ^ we see 
that the darkly shaded region is left out of the bound 
on P(decoder error |x = 00) though it is part of the 
pairwise error region. Similarly, the lightly shaded region 
is left out of the bound on P(decoder error |x = 01). 
Thus, neglecting the term P(x — » z) in the union bound 
alters the inequality. This disproves [1, Theorem 1]. The 
theorem may hold in the presence of stronger conditions 
on the constellation such as symmetry. We are unable to 
identify such sufficient conditions though. 



The problem with the proof of [1, Theorem 1] 
lies in eliminating the sub-region {r x -j n T~ z (») } where 
zW g Xc. Suppose y is received. It is reasoned in [1] 
that, if |zW — y| < |x — y|, then y should not be counted 
in the pairwise error region of (x — ► z) since the decoder 
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Fig. 2. 16QAM - The shaded regions correspond to the region not 
covered by the expurgated union bound. 



either chooses (and thus does not make an error) or 
the decoder makes an error which is already included at 
least once, in the term zW — » z. 

The right hand side of [1, (30)] can be viewed 
as an average over the transmitted signal sequence x, and 
a sum, or union bound, over z in the bad signal subset 
The decoder error region, for a given c, S and U, 
is the same for all transmitted x, and can be denoted by: 

Tc.6 = {y : min ||y - x|| 2 > min ||y - z|| 2 } 

When the union bound in [1, (30)] is expurgated, it is 
not known whether the inequality remains valid unless 
all points in r £l 6 are covered in the expurgated union 
bound for each transmitted x. So, the region r c e should 
be counted \Xg~\ = 2 d (" 1 ~ 1 ) times, each time being 
weighted by a probability measure depending on which 
was transmitted. 

Counter example to [1, Theorem 2]: Consider 
16QAM constellation with Gray labeling, as given in 
Figure |2] Let d = 1, U = and b = 0, as before. 
Let the information be transmitted in S = 2. So, some 
x £ Xq is transmitted. Suppose x = 1011 is transmitted. 
Let 2 = 0111. Setting z' = 0010 and z" = 0001, it 
can be shown that the conditions of [1, Theorem 2] 
are satisfied. It follows from the theorem that P(x — > 
z) can be neglected in [1, (30)] without altering the 
inequality. It can be similarly concluded that all z 6 Xf 
except 1111 can be neglected. Referring to Figure [2] 
it follows that the darkly shaded region is left out of 
the bound on P(decoder errorjx = 1011), though it is 
part of the pairwise error region. This occurs for all x 
of the form 10 * * (where * can be or 1). Similarly, 
the lightly shaded region is left out of the bound on 
P(decoder errorjx = 00 * *). Thus, the inequality does 
not hold and [1, Theorem 2] is disproved. 

The extension to fading channels of the above 
theorem, namely, [1, §IV.C Corollary 1], uses stronger 



conditions than the above theorem. However, it is easy to 
see how the above example works as a counter example 
to the corollary as well. 

The proof that f ex (d,n,X) is greater than or 
equal to f(d,/i,X) is hence not valid for the case of 
square QAM signal sets with Gray labeling. It is to be 
kept in mind that f ex (d, /i, X) is presented as an upper 
bound to f(d,/j.,X) in [1] only for this case - for any 
other choice of modulation, it is presented as only an 
approximation. 

III. Revised expurgated bound 

For the case of square QAM signal sets with 
Gray labeling, a revised expurgated bound, denoted by 
fex,new, is derived. Two variants of f ex ,new, denoted by 
fL,new and fex,new are presented. First, consider the 
case when d = 1. Fix a point x in the signal set, and a 
bit position i in its label. Since Gray labeling is used, the 
bit value of the corresponding position remains the same 
across either the row or the column containing x. So, if 
x £ XI , then either the entire row or the entire column 
in the signal set belongs to X£. This reduces the problem 
of identifying regions contributing to the probability of 
error to a single dimensional problem. As in any such 
single dimensional problem, the regions contributing to 
the probability of error can be covered by choosing two 
neighbors, one on each side, and constructing a union 
bound with two PEP terms. This is in contrast to the 
original expurgated bound where only one neighbor is 
considered. 

Suppose x £ XI is transmitted. Then, the 
original expurgated bound considers only the unique 
nearest neighbor in X^. An easy fix is to consider the 
nearest neighbors in X\_ h on both sides of x. This can 
be improved further by choosing the neighbors such that 
the PEP decision boundaries coincide with the actual 
decision boundaries. As an illustration, consider the 
constellation in Figure As in the counter example 
to [1, Theorem 2], suppose d = 1, U = 0, 5 = 2 
and 6 = 0. Some x £ X$ is transmitted - let 1011 
be transmitted. Since all the points in the column of 
x belong to Xq, the decision boundaries are vertical. 
The original expurgated bound includes only the unique 
nearest neighbor in Xf, namely 1111. The easy fix is to 
include the nearest neighbor in A\ 2 on the other side of 
x also - namely 0111. This variant of f e x,new shall be 
referred to as fl x new . Alternatively, the signal set can 
be extended in a lattice fashion, as shown by triangles 
in Figure and two points in the (extended) signal set 
can be identified such that the PEP decision boundaries 
coincide with the actual decision boundaries. For this 
example, the two points then are 1111 (as before) and 
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A (instead of 0111). This variant of /< 
referred to as f^ new . 

The above two methods have their relative mer- 
its and demerits. While both of the above methods yield 
upper bounds to f(l,fx,X) (this follows from a union 
bound argument), the bounds obtained using the second 
method will clearly be tighter than the first. However, 
specifying the two points is more straightforward in the 
first method. 

In either way, there are two points in X\_ h 
corresponding to each x £ X£. Let them be denoted 
by zi(x) and z 2 (x) when referring to f ex>new , by 
zf(x) and z|(x) when referring to the variant fi xnew , 
and similarly by z{ 7 (x) and z^x) when referring to 
J 'ex new Then ' fex,new can be defined by: 



shall be ip e x,new(s) given by 



fex,newiXi fti ^) 



to2 



S,U xGA-s 



In 

there 



Similar definitions hold for fl xnew and f££ new 
some cases, when a point in X£ is transmitted, 
may be no points belonging to X\_ h on one side. For 
instance, suppose some point in Xq such as 0001 is 
transmitted. Then, all points in X\ are to one side of 
the transmitted point. In such a case, zi(.) is set in the 
usual manner by choosing from these points, and Z2(.) 
is set to a special symbol H, with the understanding that 
the PEP P(x — > K) = for all x in the constellation. 
Here, H has the interpretation of a point in the extended 
constellation at infinite distance from the regular points 
of the constellation. 

The above methods are now extended to the 

g 

case a > 1 as follows. Suppose x e X^ where c = 
(ci,...,c d ), S = [h,...,id) and x = (xi, . . . , x d ). 
Each x; G X£ for 1 < I < d. Let the set consist of 
sequences of length d, where the i™ element is either 
Zi(xj) or Z2(x;). Clearly, \Z^\ = 2 d . For any two signal 
sequences, x and z, the PEP P(x — ► z) is set to zero if 
any element in z is K. Define 



fex,ne W (d, X) = }Z zZ P ^ ^ ^ 

(2) 

Since X<r is a product set, the union bound arguments 
developed for the case when d = 1 readily extend to the 
case when d > 1 to yield the following upper bound. 



f(d, fl, X) < fex,new (d, fi, X) (3) 

A computationally efficient form of f ex is derived in 
[1, (48)]. The revised expurgated bound f ex ,new can be 
expressed in a similar form, with ip e x{s) replaced by 



m2 m Z-^i L^i Z-^i 

i=l b=0 x£X? 



*A(x,»i(x))(s) + 0A(x : z 2 (x))(s)} 



The asymptotic behavior of f ex at large SNR (a « 1) 
in the presence of fading is given in [1, (62)]. A similar 
expression can be derived for f ex ,new with dh replaced 
by dh c , where d^ 2 is given by 
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|x - zi(x) 



|x - z 2 (x) 



In the design guidelines listed in [1, §V], the harmonic 
mean square distance d\ should be substituted with d\ . 
As an illustration, Table I of [1] is given here with 
the revised harmonic mean square distances. Here, d ! h 
corresponds to fi x . new , and dj/f corresponds to f^ new . 



TABLE I 

Values of d\ for some signal sets with average energy 

NORMALIZED TO 1. 
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IV. Numerical results 

In this section, the original expurgated bound 
is compared numerically with the two versions of the 
revised expurgated bound for a Rayleigh fading channel 
(K = 0) with full CSI at the receiver. The modulation 
schemes considered are 16QAM and 64QAM with Gray 
labeling. The binary code used is the standard rate- 
1/2, 64-state binary convolutional code with generators 
(ol33, ol71) used in [1] (also given in [2, pp. 507] 
as the (o634, o564) code). This code has a minimum 
distance di = 10. The revised versions of the expurgated 
bound are numerically evaluated on the same lines as the 
original expurgated bound. 

In Figure [3] the bounds on BER are graphed 
along with simulation results, for 16QAM and 64QAM 
with Gray labeling. Curves marked by 'EX orig' denote 
the original BICM expurgated bound, 'EX newl' the 
bound on BER corresponding to fl x new , 'EX new2' 
the bound corresponding to fl xnew and SIM computer 
simulation. The simulation results are obtained by using 
the suboptimal branch metric [1, (9)]. The following 
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Fig. 3. BER of BICM obtained from the optimal 64-state rate-1/2 
code and QAM signal sets. Rayleigh fading with perfect CSI. 



observations can be made from the figure. The 'EX 
newl' upper bound is greater than 'EX orig' by a factor 
of about 2 for 16QAM and about 3 for 64QAM. For 
both 16QAM and 64QAM, the 'EX new2' upper bound 
is nearly indistinguishable from (but is tighter than) the 
original expurgated bound (EX orig) for moderate to 
high SNR. This is related to djj nearly coinciding with 
dh for the square QAM constellations with Gray labeling 
listed in Table U While this suggests that the original 
expurgated bound may be a valid upper bound for square 
QAM with Gray labeling, we do not have a proof to 
support such a claim. 
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